The sensual Apollonian circle packing 
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Abstract. The curvatures of the circles in integral Apollonian circle 
packings, named for ApoUonius of Perga (262-190 BC), form an infinite 
collection of integers whose Diophantine properties have recently seen a 
surge in interest. Here, we give a new description of Apollonian circle 
packings built upon the study of the collection of bases of inspired 
by, and intimately related to, the 'sensual quadratic form' of Conway. 



1. Introduction 
In their delightful monograph entitled The Sensual Quadratic Form [8], 




The elements of P^(Q) can be viewed as the primitive vectors^ of Z^, 
considered up to sign, so that u in the image above represents ±u; such an 
equivalence class is called a lax vector. These lax vectors label an infinite 
froth of planar regions demarcated by an infinite tree of valence three. More 
precisely, this branching topograph is the graph whose edges correspond to 
bases of and whose vertices correspond to superbases, i.e. triples of lax 
vectors, any two of which form a basis. The fact that an edge connects 
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two vertices reflects the fact that any basis is contained in exactly two 
possible superbases. That is, {±-u, ±v} is contained in {±-u, ±v, ±{u + v)} 
and {±M, ±v, ±{u — v)}. One then shows that the graph is a single infinite 
valence-three tree which can be made planar in such a way that it breaks 
up the plane into regions, the boundary of each consisting of exactly those 
edges and vertices containing a single lax vector: we label the region with 
this vector. In the picture above, all the regions are labelled in terms of one 
arbitrary choice of basis, u and v. 

Just as a linear form is determined by its values on a basis, a binary 
quadratic form is determined by its values on a superbasis. The topograph 
provides lovely visual evidence of this fact for forms on Z^. Evaluate the 
quadratic form / on the lax vector labelling each region. The parallelogram 
law, 

f{u + v) + fiu-v)^2f{u) + 2f{v), 

relates the values on the four regions surrounding an edge of the topograph. 

Hence knowing the values surrounding any one vertex allows one, iteratively, 
to deduce every other value in the topograph. 

Imagining / as altitude, Conway classifies integral quadratic forms (as 
indefinite, positive definite, and so on) by their topographical terrain, de- 
scribing their lakes, wells, weirs and river valleys. 

A different, but equally beautiful picture named for geometer ApoUonius 
of Perga (262-190 BC), looks something like this. 




To form such an Apollonian circle packing, one starts with three mutually 
tangent circles of disjoint interiors (sometimes the "interior" of a circle may 
be defined to be the outside). Given such a triple, there are exactly two 
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ways to complete it to a collection of four mutually tangent circles, again 
of disjoint interiors, which is called a Descartes quadruple. In the following 
picture, a triple is drawn in solid lines, and the two possible completions 
are dotted. 




To form an Apollonian circle packing, start with any three mutually 
tangent circles, and add in the missing completions, thereby growing the 
set of circles. Repeat this process ad infinitum, at each stage adding all the 
missing completions of all triples in the collection. 




The interest for the number theorist lies in the quadratic Descartes re- 
lation between the curvatures (inverse radii) a, b, c, d of four circles in a 
Descartes quadruple: 



(For circles whose "interior is outside," curvature is negative.) This relation, 
which has been traced by Pedoe [1'^] to Descartes, in a letter to Princess 
Elisabeth of Bohemia [10, p. 49], entails that the two possible curvatures d 
and d' which complete a triple of curvatures a, 6, c to a Descartes quadruple 
satisfy the linear Descartes relation: 



In particular, an Apollonian circle packing seeded with a Descartes quadru- 
ple of integer curvatures a, b, c, d consists entirely of circles of integer cur- 
vature! Such a packing is called integral. Furthermore, it is primitive if its 
curvatures have no common factor. 

The pheromones exuded by such a mysterious collection of integers have 
inevitably attracted number theorists, with some remarkable results. See 
[3, 4, 11, 15, 19] for theorems on the density of the set in the integers, the 
occurrence of prime curvatures, and much more. 

Perhaps the reader has, at this point, noticed a similarity in the stories 
of Conway and Apollonius. In the first story, each basis (pair of vectors) 
can be completed to a superbasis (triple of vectors) in exactly two ways, 






(2) 
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Figure 1 . A primitive integral Apollonian circle packing in 
which each circle is labelled with its curvature. The outer 
circle has curvature —6. 

and the values of / on the two new possibilities sum to twice the sum on 
the original pair: 

f(u + v) + f{u-v)=2{f{u) + f{v)). 

In the second story, each triple of tangent circles can be completed to a 
Descartes quadruple in exactly two ways, and the curvatures of the two 
new possibilities sum to twice the sum on the original triple: 

d + d' ^ 2(a + b + c). 

Just as Conway labels the topograph with a quadratic form whose values 
satisfy the parallelograph law, it is natural to ask if there is a corresponding 
'Apollonian topograph,' which we can label with some type of form to 
obtain the curvatures of an Apollonian circle packing. We aim to describe 
just such an object: where Conway's topograph was essentially a picture of 
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P"'^(Q), we describe an object intimately related to P^(Q(i)). Where Conway 
labelled the topograph with values of a binary quadratic form, we consider 
Hermitian forms. 




Figure 2. A portion of an Apollonian superpacking (red), 
together with its dual (blue). The image shows those circles 
of curvature less than 20 within a square with side-length 2, 
centred on 1 + i. See Section 13. 

The 'Apollonian topograph' described in this paper is a forest. It consists 
of infinitely many trees of valence four, each of which is, in a concrete 
sense, an Apollonian circle packing: its vertices correspond to Descartes 
quadruples, and two quadruples are adjacent when they share three circles 
in common (Theorem 4.1). The Apollonian forest, in a different form called 
the 'Apollonian city' (formed of infinitely many 'palaces'), can be labelled 
with the imaginary part of any Hermitian form. For a particular choice of 
form, we obtain the curvatures of the circles in the corresponding packings; 
it is also possible to obtain the centres in this manner (Theorem 4.2). 

With this perspective comes a natural way to associate a full rank sub- 
lattice of T? to each circle in any Apollonian circle packing, and there is 
a corresponding 'lattice Descartes rule' (Theorem 4.5) which describes the 
lattices of the two completing circles for a triple. To each circle we also 
discover an associated equivalence class of quadratic forms, whose values, 
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translated, give the curvatures of its tangent circles (Theorem 4.6). That 
such forms existed was known [3, 11]; here we obtain a precise description 
of these forms in terms of circle lattices, which, in particular, describes the 
form of a circle in terms of the forms of its neighbours (Theorem 4.7). 

Taking the entire forest as one object, we obtain a union of packings 
called an Apollonian superpacking, first introduced in [ ], and shown in 
Figure 2. The superpacking contains copies of every possible primitive 
Apollonian circle packing. Through its connection to the Apollonian forest, 
the superpacking can be viewed as a beautiful picture of the collection 
of full-rank sublattices of Z^, a picture of PGL2(Z[?])/ PSL2(Zi), or, from 
another perspective, a piece of a picture of P^(Q(i)). Another piece of the 
picture of P^(Q(i)) inspired by Conway's topograph was studied by Bestvina 
and Savin in [2] (see Section 17). 

Acknowledgements. The author would like to thank Lionel Levine for draw- 
ing her attention to Apollonian circle packings, and would especially like 
to thank David Wilson for sharing data he had collected associating cir- 
cles and lattices in the context of abelian sandpiles. It was an examination 
of his data that led to this study. For more on this connection, see [17]. 
The author owes a debt of gratitude to the invaluable and detailed sugges- 
tions provided by Jeffrey Lagarias, Andrew Granville, Jonathan Wise and 
an anonymous referee on an earlier draft. The images in this paper were 
produced using IPE [ ] and Sage Mathematics Software [20]. 



2. The Apollonian city 

Our central object is the Apollonian city and its labelling by lax vertex 
lattices. After defining these in two brief sections, we can then state the 
Main Bijection relating the Apollonian city to Apollonian circle packings. 
Later, we will see how this labelled object is the analogue of Conway's 
topograph. 

Extending terminology from Conway and Fung's monograph, a lax vector 
is an equivalence class {u, —u,iu, —iu} for some primitive^ u G Z[z]^. A 
superhasis'^ is a collection of three lax vectors {u.,v,w}, any two of which 
form a basis for Z[i]^. This entails that eiu + + e^w = 0, for some 
appropriate choice of units e^. 

The 2-cells of the Apollonian city will consist of the collection of all 
superbases of Z[i]^, in the form of triangles whose sides are labelled with 
the superbasis elements u, v and w: 



^i.e. not a Gaussian integer multiple of another vector of Z[i]^. 

''if we consider P-'^(Q(i)) as P^(Z[i]) over SpecZ[i], then a superbasis can be defined 
as three distinct Z[z]-points of P^(Z[z]). 
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W 



We will call the 2-cells walls. 

An ultrabasis is a collection of six lax vectors {u,v,w,x,y, z} such that 
the four subsets {u,v,w}, {u,y,z}, {v,x,z} and {w,x,y} are superbases 
(as it turns out, if a set of six lax vectors is an ultrabasis, it is an ultrabasis 
in exactly two ways; we consider these distinct). Whenever the union of four 
superbases forms an ultrabasis, the sides of the corresponding triangular 2- 
cells are identified according to their labels, and the Apollonian city includes 
a tetrahedral 3-cell whose boundary consists of the four triangles: 




W 



We will call the 3-cells chambers. 

This information is sufficient to define the Apollonian city as a 3-dimensional 
cell complex whose edges are labelled with lax vectors: simply glue the col- 
lection of all 2-cells according to their formation into superbasis 3-cells. 
According to the construction, a given lax vector may appear as a label 
on more than one edge, but a given superbasis or ultrabasis has a unique 
appearance as a 2-cell or 3-cell respectively. 

Each component of the Apollonian city is called a palace. 

We will see that each chamber is adjacent to exactly four others, sharing 
one wall with each; the valence-four graph whose vertices are chambers 
connected according to adjacency is called the Apollonian Forest (described 
in more detail in Section 5); we will see that it consists of infinitely many 
trees. 
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3. The Vertices of the Apollonian city 

Suppose we have a superbasis represented by vectors u,v,w with the 
convention that u + v + w = (i.e., representatives u, v and w are chosen so 
that no unit multiples are required to cause the sum to vanish). Then there 
are exactly two ways to form an ultrabasis containing this superbasis'': 




w w 
For example, the leftmost triangle is a superbasis because 



u — i{u + iv) + {w + iu) = u + v + w = Q. 

These two ultrabases share the face u,v,w in the Apollonian city. Any 
ultrabasis is of the form of one of these two standard ultrabases, for some 
u, V, w; and any pair of chambers sharing a wall in the city is of the form of 
this pair. 

We will now label each vertex of the standard ultrabases above with an 
ordered pair of vectors, according to the following picture: 



{u,iv) {u,iv) 




{w,iu) {v,iw) {w,iu) {v,iw) 

Definition 3.1. Considering as a Z-module of rank four, let L be an 
oriented Z-submodule of rank two generated by an oriented Z[z]-basis for 
Z[z]^. The lattice L is called a vertex lattice and the pair {L,iL} is called 
a lax vertex lattice. 



Note that the three edges touching a single vertex in an ultrabasis also form a su- 
perbasis (see Section 17 for more about this 'duality'). 
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For example, the lattice {a(— z, — 2 — z) + b{—l,—l + i) : a, 6 G Z} is 

f -i -1 \ 

a vertex lattice, since the matrix „ . ^ . has determinant — 1. 

\^— 2 — I —1 + iJ 

This vertex lattice, together with the vertex lattice {a(l, 1 — 2i)+6(— z, — 1 — 
i) : a, 6 G Z} form a lax vertex lattice. 

We will attach to each vertex of the Apollonian city the lax vertex lattice 
generated by the labels in the picture above (i.e. label each ultrabasis as a 
standard ultrabasis for some u, f, w, and then attach labels as shown). 

Proposition 3.2. The lax vertex lattice attached to a vertex in the Apollo- 
nian city is well-defined. 

The proof can be found in Section 7. 



4. From a city to a packing: main results 

The space P^(C) is defined as the collection of non-zero vectors in 
modulo multiplication by non-zero scalars. We identify P^(C) with points 
of C U {oo} by the map (a, h) i— )■ a/h. The subset P^(Q(z)) consists of those 
equivalence classes of vectors containing a vector in Q(z)^; one may identify 
this with the collection of primitive vectors of Z[i]^ up to sign (i.e. lax 
vectors) . 

Mobius transformations act on P-'^(C) via the matrix action of GL2(C) 
on vectors in C^. Since matrices which are scalar multiples have the same 
action, we can may consider instead the action of equivalence classes of ma- 
trices up to non-zero scalar multiple; this is the group PGL2(C). The group 
PGL2(Z[i]) is taken to be the subgroup of PGL2(C) formed of equivalence 
classes containing a matrix with entries in Z[i] and determinant a unit in 
Z[z]. Those with such a representative whose determinant is ±1 form the 
smaller subgroup PSL2(Z[i]). Finally, within this group, some may have a 
representative with entries in Z and determinant 1; this is PSL2(Z). We 
obtain a chain of subgroups: 

PGL2(C) > PGL2(Z[?]) > PSL2(Z[i]) > PSL2(Z). 

In what follows we will frequently refer to an equivalence class of PGL2(C) 
by a matrix representative. The group PGL2(Z[i]) takes elements of P"'^(Q(i)) 
to itself, so it can be considered to act on the lax primitive vectors of Z[z]^, 
which will sometimes be identified with elements of Q(i) U {oo} via the map 
(a, h) H- >■ a/6 above. 

As described in the introduction, an Apollonian circle packing is a col- 
lection of infinitely many circles in P^(C) = C U {oo} (where lines are 
considered circles passing through oo) with disjoint interiors, formed from 
an initial triple of three mutually tangent circles by an iterative process. 
At each stage, for every triple of mutually tangent circles already included, 
one adds to the collection the two circles which complete that triple to a 
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Figure 3. An Apollonian circle packing (in black) with the 
Apollonian city (in grey) shown as an overlay; the vertices 
of the city are located at the centre of the corresponding cir- 
cles. (Only circles and vertices of the Apollonian city with 
curvatures below a small bound are shown; the vertex corre- 
sponding to the outermost circle is ommitted.) 

Descartes quadruple (we will call these the completing circles). In the limit, 
we obtain a well-defined infinite collection of circles. 

Such a packing, or a subset of one, is called integral if all of its curvatures 
are integral, and strongly integral if, in addition, the centre of each circle, 
multiplied by the curvature, is a Gaussian integer. It is called primitive if 
the collection of curvatures has no non-trivial common factor. 

Each circle comes with a designated interior; equivalently, we can consider 
the circle oriented so that the interior is the region to the right as one travels 
around the circle in the direction of orientation. 

The central theorem (the 'Main Bijection') describes the relationship be- 
tween the Apollonian city and Apollonian circle packings: we show that the 
Apollonian city comes in infinitely many connected components, or palaces, 
the vertices of a particular palace being in bijection with the circles in a 
particular primitive strongly integral Apollonian circle packing, where ad- 
jacent vertices correspond to tangent circles. Placing each vertex of the city 
at the centre of the corresponding circle of the packing, we obtain pictures 
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like Figure 3. As a consequence of the Main Bijection, the Apollonian city 
is an example of a class of graphs named Apollonian networks after this 
connection [i]. 



Theorem 4.1 (Summary of Theorems 8.1 and 12.1). We have aPGL2{1j[i])- 
equivariant bijection (the 'Main Bijection') between the following sets: 



Under the bijection, circles which are tangent within an Apollonian circle 
packing correspond to vertices connected by an edge in the Apollonian city, 
and these properties are preserved by the action. In particular, the vertices 
of a chamber in the city correspond to a primitive strongly integral Descartes 
quadruple of circles. 

Furthermore, the Apollonian city consists of infinitely many components, 
each of which corresponds to a unique primitive strongly integral Apollonian 
circle packing. Finally, every primitive Apollonian circle packing appears, 
up to rigid motions. 



The bijection is mediated through the lax vertex lattice: its basis vectors 
form the columns of the corresponding element of PGL2(Z[i]). 

It may be informative to carry along an example of this bijection. Choos- 
ing u = {i, 1) and v = (0, 1), by convention we set w = {—i, —2) (so that 
u + v + w = 0). We obtain the following ultrabasis tetrahedron, with vertex 
lattices shown (as matrices whose column vectors generate the lattice): 



< 



oriented circles and lines in C 
which are images of the real line 
under Mobius transformations 



> 




12 



KATHERINE E. STANCE 




The corresponding images of M in Si form a Descartes quadruple: 






The circles of Si fill the complex plane densely. Figure 2 shows those 
circles that have curvature less than 20. The red circles, no two of which 
cross transversely, correspond to the subgroup PSL2(Z[z]) of PGL2(Zi[i]); 
the blue circles correspond to its nontrivial coset. 
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Just as Conway's topograph could be labelled with various quadratic 
forms, we will label the Apollonian city with the imaginary part of a Her- 
mitian form at each vertex (evaluated via the vertex lattice). These values 
satisfy the linear Descartes rule, and for a special choice of Hermitian form, 
we obtain the curvatures, and the products of the centres and curvatures 
of the circles. 

Theorem 4.2 (Summary of Propositions 10.2 and 11.1). Let H : CxC — > M 
be the imaginary part of a Hermitian form. Then H takes a well-defined 
value on a lax vertex lattice, and hence on a vertex of the Apollonian city, 
where this value is obtained by evaluating H on any Z-basis of the lattice. 
Consider the values A, B and C of H on the vertices of the common wall 
between a pair of adjacent chambers in the Apollonian city. Let D and D' 
be the values at the other two vertices of the chambers. Then 



In particular, consider the following imaginary parts of Hermitian forms: 



The value of Hi is the curvature of the corresponding circle under the Main 
Bijection, and H2 + iH^ is the curvature times the center of the circle. 

Continuing the example quadruple above, we find the centres and curva- 
tures of the Descartes quadruple as follows: 



D + D' = 2{A + B + C). 



Hi{{a,f3),{j,6)) 
H2{ia,P),i^,5)) 
H,{{a,P),{^,6)) 




.5), 



curvature = -1 
center = i 




1 + i 
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Because Hi of Theorem 4.2, whose values give curvatures, depends only 
on the second coordinates of its vector inputs, one may choose to restrict 
attention at a vertex of the Apollonian city from its vertex lattice to the 
Z-module generated just by the second coordinates of the vertex lattice. 

Definition 4.3. Considering Z[i] as a Z- module of rank two, let L be a 
full-rank Z-submodule. The lattice L is called a circle lattice. If L has basis 
(a, b) and (c, d), then the pair {L, L'} where L' has basis (6, —a), (rf, — c), is 
called a lax circle lattice. 

For example, the vertex lattice {a{—i, — 2 — i) + 6(— 1, —l + i):a,bE Z}, 
corresponds to the circle lattice {a(— 2, —1) + 6(— 1, l);a,b G Z}. 

Thus we associate to each vertex of the Apollonian city a lax circle lattice. 
This is no longer a bijection; each lax circle lattice may appear many times. 

Proposition 4.4 (See Proposition 15.1). The determinant of a circle lattice 
is the negative of the curvature of the associated circle. 

Just as there is a Descartes rule for curvatures, we can now record a 
'Descartes rule for circle lattices'. 

Theorem 4.5 (Theorem 15.2). Consider three mutually tangent circles in 
an Apollonian circle packing. For some integers a, b, c, d, they have the 
following circle lattices (represented by matrices whose column vectors form 



a basis): 




The two completing circles have the following circle lattices: 




Figure 4 shows the circle lattices associated to the circles in an example 
Apollonian circle packing. In the example we have been following, the circle 
lattices are as follows: 
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The relationship to Conway's topograph and the Apollonian city extends 
beyond an analogy: one obtains infinitely many copies of the topograph (in 
a slightly different form) within the Apollonian city. 

Theorem 4.6 (Introductory form of Theorem 14.1). Fix a vertex of the 
Apollonian city with vertex lattice L. Then the adjacent vertices may he put 
in bijection with the primitive elements x & L in such a way that any two 
such vertices are connected by an edge of the Apollonian city if and only if 
their corresponding x,x' & L form a basis for L. Furthermore, for H the 
imaginary part of a Hermitian form, the values of H at the vertices are the 
values of a translated binary quadratic form on L. 

If H is the curvature form Hi of Theorem 4.2, then the quadratic form 
described in this theorem is given by the circle lattice of the fixed vertex 
in a simple way: if M is a matrix representation of the circle lattice, then 
M^M is the matrix of the quadratic form. As a consequence, the form has 
discriminant negative four times the square of the curvature of the circle. 

Theorem 4.7 (Corollary 16.2). We have an explicitly computable bijection: 



To any circle in an Apollonian circle packing are associated a lax circle 
lattice and hence a quadratic form; the values of the form, translated by 
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Figure 4. An Apollonian circle packing with circle lattices. 
(Only circles with curvatures below a small bound are shown.) 

the determinant of the lax circle lattice, are exactly the curvatures of all 
the circles in the Apollonian circle packing which are tangent to the given 
circle. 

For example, fixing the lower circle of curvature 2 in Figure 4 gives rise 
to the quadratic form + AY'^ . The curvatures of the tangent circles are: 

-1,2,3,6,11,15,18,... 

which list, when translated by 2 (the curvature of the fixed circle), are 
exactly the values 

1,4,5,8,13,17,20,... 

of the quadratic form. 

This recovers a relationship between Apollonian circle packings and qua- 
dratic forms that has already been studied; see for example [12]. In fact, 
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it is one of the main ingredients in the proof of several major Diophan- 
tine results concerning the curvatures of an Apollonian circle packing [•!, 4]. 
Through Theorems 4.5 and 4.7, we obtain a simple new rule for relating 
the quadratic forms of tangent circles. 



5. The Apollonian forest 

The Apollonian forest, shortly to be defined as a type of dual to the 
Apollonian city, is at times more convenient, and it is more clearly an 
analogue to Conway's topograph. 

The forest is a graph whose vertices and edges are the collection of ul- 
trabases and superbases of respectively. Their adjacency is defined 

by the property of an ultrabasis containing a superbasis as one of its four 
defining superbases. Hence we obtain a valence-four graph. A small piece 
of the Apollonian forest looks like this, where (U) represents an ultrabasis 
vertex, and — ^ represents a superbasis edge: 



■<u>- 



■<u>- 



■<u>- 



■<u>- 



Superimposing the forest on the city, the forest's vertices can be consid- 
ered to lie in the center of the corresponding chambers of the city, and its 
edges to pierce the walls of the chambers (this is the 'duality' referred to 
above) . 

Once we have established the Main Bijection, each component of the 
Apollonian forest will be in bijection with an Apollonian circle packing; 
the tree is in fact the graph whose vertices are the Descartes quadruples of 
the packing, and whose edges connect pairs of Descartes quadruples which 
share a common triple of circles. 

The Apollonian forest is the true analogue of Conway's topograph. Sup- 
pose we form a 2-dimensional cell complex based on the topograph, where 
each face gives rise to a vertex, and faces which touch become vertices 
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connected by an edge. Each vertex of the topograph gives rise to a 2- 
dimensional triangular face. In this new picture, each vertex is labelled 
with a vector up to sign, each edge denotes a basis, and each face (triangle) 
denotes a superbasis. This new object, shown in black against the grey 
topograph in the following picture, is the analogue to the Apollonian city. 




It is possible, although we do not prove it here, to embed the Apollonian 
forest in three-dimensional space in such a way that it breaks the space up 
into regions, each corresponding to a lax vertex lattice, i.e. a vertex of the 
Apollonian city. 

6. The Apollonian forest consists of trees 

The following will be needed to prove Proposition 3.2, among others. 

Proposition 6.1. Every component of the Apollonian forest is a tree. 

Proof. Consider the form / : Z[i]^ — )■ Z given by 

f{a + bi,c + di) = + 1?. 

This induces a labelling on the vertices (ultrabases) of the forest as follows: 
label an ultrabasis of the form 




w 

with the integer 

f{u) + f{v) + /H + f{x) + f{y) + f{z), 
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which we will call the height of the ultrabasis. This is well defined since 
f{u) = f{-u) = f{iu) = f{-iu). 

We will shortly show that at any vertex, at least three of the four adjacent 
vertices have a strictly greater height. If this is the case, then, starting 
out at any vertex and travelling upward, we must travel upward at every 
subsequent step, unless we double back. From this, one concludes that there 
are no non-trivial loops. 

We consider one of the standard ultrabases and its four neighbours (in 
particular, u + v + w = 0), shown in Figure 5. 




w 



Figure 5. A portion of an Apollonian tree. 

Suppose that u has first coordinate a+bi, v has first coordinate c+di, and 
w has first coordinate {—a — c) + {—b — d)i. Then the height of the central 
ultrabasis minus the heights, respectively, of each of the four neighbours is, 
clockwise from top: 



- 12{ad-bc), 

- 12((a + + {b + df -ad + be), 



-12(c^ + c/^ - ad + bc), 
-12{a^ + b^ - ad + bc). 
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Since u,v form a basis for Z[z], their first coordinates, a + bi and c + di, 
must be coprime, hence ad — bc^ 0. If the first of the four quantities hsted 
above is non-negative, then all the others are negative. In other words, if 
one of the four neighbours has an equal or lower height than the central 
ultrabasis, then the other three are higher. The proof for the other standard 
ultrabasis is much the same. Without loss of generality, this proves the 
general case. □ 

We will see later that the Apollonian forest consists of infinitely many 
trees. 



This section contains the proof of Proposition 3.2 showing that the vertex 
lattices associated to a vertex of the Apollonian city are well-defined; as 
a consequence of the proof, we also discover that each lax vertex lattice 
appears exactly once in the Apollonian city. 

Proof of Proposition 3. 2. As a means of proof, we will give an equivalent 
definition of the labelling of the vertices which is clearly well-defined. This 
involves adding a little extra structure to the labelling of the Apollonian 
city. 

Given an edge in the city labelled with u, representing an equivalence 
class of vectors modulo units, we instead label each end of the edge (i.e. 
we label the edge together with a choice of one of its two vertices; we 
place the label closer to the chosen vertex) with an equivalence class of 
vectors modulo sign. The two sign-equivalence classes have as their union 
the original unit-equivalence class, so that the edge 



There are, of course, two ways to do this (differing by multiplication by i). 

This splitting of labels can be performed for every edge of a tetrahedral 
chamber in such a way that surrounding a vertex, the three nearby labels 
X, y and z satisfy x + y + z = for some choice of signs (instead of unit 
multiples). Up to multiplication of every label by i, there is a unique way 
to do this. For example, here is a standard ultrabasis tetrahedron before 
and after splitting (signs are suppressed). 



7. Lax Vertex Lattices 




becomes 
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w w iw 



Since the forest has no loops (Proposition 6.1), this shows that, once 
a single label has been split, there is a unique way to split labels on all 
the ultrabases in a single palace simultaneously and consistently. Multi- 
plying every label simultaneously by i allows us to move from one of the 
two possible global splittings to the other. (Choosing one of these two 
global splittings is an arbitrary choice we will make for the moment; this 
discomfiture will be resolved shortly.) 

Now each vertex in the Apollonian city has many proximal labels on its 
adjacent edges, but any three arising from one tetrahedron satisfy x+y+z = 
for some choice of signs. This implies that, having selected two such 
proximal labels (from a single tetrahedron among those sharing that vertex), 
say u and iv, all other proximal labels are Z-linear combinations of u and iv. 
This defines a rank two Z-submodule of where the latter is considered 
as a free Z-module of rank 4; this will be the vertex lattice. The submodule 
is generated by two primitive vectors forming a Z[z]-basis for 

All that remains is to explain how to determine the orientation of the 
vertex lattice. Any triangular side which touches the vertex has proximal 
labels X and y and distal labels ix and iy (up to sign) on the two edges 
connected to that vertex. Choose the signs of the labels on the entire 
triangle so that two things hold: 

(1) if multiplication by —i is considered to orient an edge, the triangle 
is a cycle, 

(2) the 'forward' labels on the three edges sum to zero. 

For example, if u + v + w = 0, then there are exactly four ways to do this: 




w iw —w —iw —w iw w —iw 



At this point the labels represent vectors, not vectors up to unit or sign. 
Now, put an order on the basis generating the submodule of our chosen 
vertex according to the orientation of the triangle. In this case, the ordered 
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basis at the upper vertex is {u,iv), {—u, —iv), {iv, —u) or {—iv,u), for the 
four pictures, respectively, but they all represent the same oriented lattice. 
It remains to verify that if we do this with another wall of the same chamber 
(under the same global splitting), we obtain the same oriented lattice at a 
common vertex. It suffices to check this for one other wall of a standard 
ultrabasis^: 




w + iu iw — u 



At this point we have now defined a vertex lattice attached to the specified 
vertex of the Apollonian city. However, we made one arbitrary binary choice 
in the process, resulting in two possible global labellings of the vertices of 
each palace, differing by multiplication by i. 

Hence, what we have actually defined, without making an arbitrary 
choice, is a lax vertex lattice. 

Finally, the lax vertex lattices associated to the standard ultrabases are 
exactly those defined in Section 3. □ 

In some circumstances, for a single palace, it will be convenient to make 
the arbitrary binary choice described in the proof above, and talk about 
the (non-lax) vertex lattice at a vertex. 

We can summarize the 'meaning' of the various pieces of the Apollonian 
city or forest in a table. 



Apollonian forest 


Apollonian city 


meaning (up to unit) 




vertex 


lax vertex lattice 




edge 


primitive vector of Z[i]^ 


edge 


triangle (wall) 


superbasis of Z[i]^ 


vertex 


tetrahedron (chamber) 


ultrabasis of Z[i]^ 



Proposition 7.1. With the lax vertex lattice labelling described above, each 
lax vertex lattice appears exactly once in the Apollonian city. 

Proof. For every basis m, f of Z[2]^, there are exactly four possible lax vertex 
lattices (exemplified by bases (m, f), {u,iv), {u, —v) and (m, —iv)). All four 
appear somewhere in the Apollonian city. For example, the lattice generated 
by M, iv appears on one vertex of the superbasis m, f , —u—v, while the lattice 
generated by u, v appears on one vertex of the superbasis u, —iv, —u + iv. In 
fact, each lax vertex lattice appears only once, since the Z-bases {u,v) are 

^Notc that this sign-refinement of the spht labelhng and orienting of edges cannot be 
performed in a single consistent way simultaneously throughout a palace. But it is not 
necessary to do so. 
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Si 



in bijection with the superbasis walls u, iv, —u — iv which share that vertex 
(this follows from the proof of Proposition 3.2), and hence the lax vertex 
lattice uniquely determines the collection of superbasis walls it touches. □ 

8. The main bijection 
Theorem 8.1. Consider the following sets with actions by PGL2(Z[i]).- 

oriented circles and lines in C 
which are images of the real line 
under Mobius transformations 

^ M /or M G PGL2(Z[i]) 

with action by left multiplication on M . 

(2) S2 = PGL2(Z[?])/PSL2(Z) with action by left multiplication. 

(3) S3 = {lax vertex lattices} with action on the lax vertex lattice. 

(4) = {vertices of the Apollonian city} with action on the lax vertex 
lattice labelling the vertex. 

These four sets with action are in PGL2{'Ij[i])-equivariant bijection. In par- 
ticular, M E S2 corresponds to 

(1) the image of M in Si, 

(2) the lax vertex lattice generated by the column vectors of M in S3, 

(3) the vertex labelled with the lax vertex lattice generated by the column 
vectors of M in S^. 

Proof. Represent a lax vertex lattice by a matrix M whose column vectors 
form an ordered Z-basis for the lattice. Matrices whose column vectors 
form a basis for the same lax vertex lattice must be considered equivalent. 
In particular, 

(1) M is equivalent to MN for N e PSL2{Z), 

(2) M is equivalent to iM, -M and -iM. 

Hence, a lax vertex lattice is in bijection with an element of 

PGL2(Z[z])/PSL2(Z) 

(here the quotient is by right multiplication). 
The group of Mobius transformations PGL2(Z[z]) takes circles to circles in 
(C) (where as usual, straight lines are considered circles passing through 
00). Since PSL2(Z) is the group of Mobius transformations in PGL2(Z[z]) 
that fix the real line, preserving its orientation, to specify an element M G 
PGL2(Z[i])/ PSL2(Z) is to specify an oriented circle which is the image of 
the real line (oriented to the right, say) under M. 

Clearly, these bijections preserve the actions of PGL2(Z[i]). Finally, we 
have already seen that lax vertex lattices are in bijection with vertices of the 
Apollonian city (Proposition 7.1). We have defined the action of PGL2(Z[2]) 
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on the vertices by its action on the vertex lattices, so the action is preserved 
by the bijection. □ 

9. The action of PGL2(Z[z]) on the Apollonian city 

Proposition 9.1. The action o/PGL2(Z[i]) in Theorem 8.1 preserves ad- 
jacencies in the Apollonian city, or equivalently, the Apollonian forest. The 
Apollonian forest consists of infinitely many trees, and the subgroup taking 
any one tree to itself is a conjugate of 

^-((: :l))<p« 

Proof. The group PGL2(Z[z]) acts on P^(Q(i)) in such a way that a super- 
basis is taken to another superbasis, and an uhrabasis is taken to another 
ultrabasis. This gives an action of PGL2(Z[i]) on the Apollonian forest, 
or equivalently, on the Apollonian city, preserving adjacencies. This action 
is conjugate to the action in Theorem 8.1 on the vertices of the Apollo- 
nian city, because the action of M on a basis u, v corresponds to the action 

of ( „ . M \ ^ . 1 on the lax vertex lattice with ordered basis u, iv. 



ij \0 i^ 

The remainder of the proof will discuss this action in lieu of the action in 
Theorem 8.1. It will suffice that the actions are conjugate. 

Consider a superbasis u + v + w = 0. The subgroup fixing this edge of 
the forest, in terms of basis u,v, is generated by the two elements 

OA p _fo -1 
1 oj' ^'-[i -1 

The matrix is of order i, and this group of order 6 corresponds to the 
permutations of the elements u, v and w. Although these elements fix the 
edge, they do not fix the entire forest. Instead, they correspond to flipping 
the edge (swapping its two vertices), and rotating the edge (acting cyclically 
on the three adjacent edges at each of its vertices). 

One edge has six adjacent edges. Representative elements moving the 
edge to each of these adjacent edges are 



Therefore, the tree in the forest containing a given superbasis consists 
of the orbit of that superbasis under the group F generated by all these 
matrices. Since F2 = F^MF^M-^F^MFs, we have 

T = {M,Fs). 
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The first generator is parabolic and the second elliptic of order 3. This 
subgroup has infinite index in PGL2(Zi[z]) Table 4], hence the Apollonian 
forest has infinitely many connected components. By Proposition 6.1, they 
are all trees. □ 

The action of a certain group called the Apollonian group on the col- 
lection of ordered, oriented Descartes quadruples of curvatures has been a 
principle tool in the study of Apollonian circle packings [12, 14]. The Main 
Bijection shows that each tree in the Apollonian forest is associated to an 
Apollonian circle packing in such a way that the vertices correspond to 
Descartes quadruples and edges indicate quadruples which share a common 
triple of circles. The action described above is different than the action of 
the Apollonian group. For example, the quadruples are unordered, and the 
action above does not take every vertex to an adjacent vertex. 



10. Labelling the Apollonian city with a Hermitian form 

In order to establish the rest of the properties of the bijection in Theorem 
4.1, we must establish several results relating to Hermitian forms. We will 
describe a method to display the values of a Hermitian form on both the 
edges and vertices of the Apollonian city, in analogy to the values of a 
quadratic form displayed on Conway's topograph. 

A Hermitian form is a pairing 

/i : C" X C" ^ C 

satisfying for any u,v,w G C" and a G C, 

(1) h{au + v,w) = ah{u, w) + h{v, w) 

(2) h{u, av + w) = ah{u, v) + h{u, w) 

(3) h{u, v) = h{v, u) 

We will study the imaginary part of h (i.e. the coefficient of i). Then 
H has the following important property: 

H{au + 6f , cu + dv) = {ad — bc)H{u, v) for a, b,c,d E M. (HI) 

(Note: since ih{u, v) = h{u, —iv), we may have equivalently chosen to study 
the collection of real parts of h; in fact, h is determined by H.) 

From now on, we will assume h is integer valued on x i.e. it 

takes values in 

We can evaluate if at a vertex of the Apollonian city by applying H to 
a basis of the lax vertex lattice. Properties (HI) and conjugate linearity 
together guarantee that the if-value is well-defined. We can also evaluate 
H on an edge of the Apollonian city by applying the map H{ix, x) to the 
edge label x (also well-defined by conjugate hnearity). 

As a demonstration, here is the usual superbasis triangle labelled with u, 
V and w (where, as always, u + v + w = 0). A split labelling and orientation 
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(as in the proof of Proposition 3.2) are also shown, to demonstrate that 
they ahgn nicely with the if-labelling. 




H{iv, v) 



w H{iw,w) iw H{v,iw) 



Figure 6 shows the Hermitian labels on the first and second standard 
ultrabasis tetrahedra. 

Lemma 10.1. Ifu + v + w = Q, then the following quantities are real: 

(1) h{u, iv) + h{v, iw) — h{iv, v) 

(2) h{w — iv, iu — v) + h{w + iv, iu + v) — h{iv, v) — h{iw, w) — h{iu, u) 

Proof. The proof is by direct calculation. It is convenient to note that for 
a Hermitian form, the following quantities are real: 

h{u, u), h{u, v) + h{v, u), i{h{u, v) — h{v, u)). 

Note also that h{ix,iy) = h{x,y) and that h{ix,y) = —h{x,iy) for all 



x,y. 



□ 



As a consequence, we discover that the if-labels on the pair of standard 
ultrabase tetrahedra have the following relationships: 

Proposition 10.2. With the H-labelling described above, the following 
hold: 

• the label on an edge is the sum of the labels at its endpoints, i. e. 

H{u,iv) + H{v,iw) = H{iv,v) (H2) 




A + B = C 

• if the two standard ultrabases are considered together, the labels on 
their common edges add up to the sum of the labels of their non- 
common vertices, i.e. 

H{w — iv,iu — v)+H{w + iv,iu + v) = H{iv,v)+H{iw,w) + H{iu,u) (H3) 
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H{u, iv) 




H{w,iu) H(iw,w) H{v,iw) 

H(u, iv) 




H{w,iu) H{iw,w) H(v,iw) 

Figure 6. The first and second standard ultrabasis tetrahe- 
dra, labelled with the imaginary part of a Hermitian form. 
Note that the labels on the second are the same as the first, 
except that sums become differences and differences become 
sums. 
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D + D'^A + B + C 

if the two standard ultrahases are considered together, the labels on 
their common vertices, doubled, add up to the sum of the labels of 
their non-common vertices, i.e. 



H{w—iv, iu—v)+H{w+iv, iu+v) 

A 



2{H{u,iv)+H{v,iw)+H{w,iu)) (H4) 
A 




.- C 



D + D' = 2{A + B + C) 

By Proposition 7.1, the collection of Hermitian values 

{H{u,v) : u,v form a Z[z]-basis for ^[i]^} 

is now displayed on the vertices of the Apollonian city. Property (H4) can 
be considered the linear Descartes rule for Hermitian forms. 

11. The curvature and centre of a circle in P^(C) 

By oriented curvature of an oriented circle, we mean a quantity that is 
positive if the circle is oriented clockwise, negative if oriented counterclock- 
wise, and is equal to the usual curvature in absolute value. 



Proposition 11.1. Let 



M 



V2 






1 

V2 



a 7 
P 6 



Then the oriented curvature of the oriented circle which is the image of the 
real axis (oriented to the right) under M , is given by 

Im {/36) . 
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Furthermore, the product of the curvature and centre of the circle is given 
by 

i{^(3 — a5) = Im (/37 + a6) + ilm (i(3^ — ia6) 

Furthermore, this value has even real part and odd imaginary part whenever 
M G PSL2(Z[i]); it has odd real part and even imaginary part whenever 
M G PGL2(Z[z]) \ PSL2(Z[2]). 

Finally, the circles in Si of Theorem 8. 1 have integer curvature, and the 
product of curvature and centre is a Gaussian integer. 

Proof of Proposition 11.1. Define functions on elements = ( ^ ] ) of PGL2(Z[«]) 
by 

Hi{N) = Im(2/35), 

H2{N) = lm + , 

HsiN) = Im {i(3^-ia6) . 

Then each is the imaginary part of a Hermitian form on the column 
vectors of the matrix (note that Hi is twice the form in the statement 
of the Proposition). We will show that the image of the real line under N 
has curvature Hi{N) and center x curvature H2{N) +iH3{N), from which 
the statement of the Proposition follows (since the image of M is the circle 
obtained from the image of by dilating the complex plane by a factor of 
two). 

Each Hi, and the image circle itself, is invariant under replacing N by 
any other transformation N' = NP for some P G PSL2(M) (see property 
(HI) of the imaginary part of a Hermitian form in Section 10). If we change 

the orientation of the circle, replacing with A^ ^ ^''^ , then the values 

of the Hi change sign. 

Keeping track of the possible change of orientation and sign, we can 
replace A^ by a transformation of the following form (using an appropriate 

P G PSL2(M) and possibly ^ ^ ^ as above), whose image is the same 

circle: 

' k+b 

j^'=lV2k "v^ 1 G PGL2(C), beQ{i),ke 




2k 

The image of the real axis under this A^' is the clockwise circle with centre 
b and radius k. To verify these statements, compose the transformation 

1/v^ '''^'Jv^ taking M (oriented to the right) to the clockwise unit 
elements of PGL2(C). 
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For this A^', Hi{N') = 1/k, and b/k = H2{N') + iH^{N'). Note that 
iipfP — a5) = Im (/37 + a6) + ilm (if3^ — iaS) . 

Finally, 

a5 - 7/3 - {a6 - 7^) = a{6 -6)--f{l3-^)=0 (mod 2) 

Since a6 — 7/3 is a unit, this means the real and imaginary parts of the 
centre-times-curvature i{'yf3 — aS) consist of one even, and one odd, integer. 
In particular, if is in PSL2(Zi[i]), then the real part is even, and the 
imaginary part is odd. If is in the nontrivial coset of PSL2(Z[i]) in 
PGL2(Z[i]), then the opposite is true. □ 

12. Descartes Quadruples correspond to Apollonian 

chambers 

The goal of this section is to prove the following theorem, which will 
complete the proof of Theorem 4.1. 

Theorem 12.1. Under the bijection of Theorem 8.1, primitive Descartes 
configurations of circles in Si are in bijection with chambers of the Apollo- 
nian city. Furthermore, the palaces of the Apollonian city are in bijection 
with the collection of primitive Apollonian circle packings in Si. Finally, up 
to rigid motions, every primitive Descartes configuration of circles appears 
in Si, as does every primitive Apollonian circle packing. 

We require a sequence of propositions and lemmas. 

Proposition 12.2. The circles of Si in bijection with two adjacent vertices 
of the Apollonian city ( as in Theorem 8.1) are tangent. 

Proof. The lax vertex lattices at two adjacent vertices can be given in the 
form 

{u,iv) and {iv,u + v). 
Therefore, the corresponding circles both contain the point iv (considering 
a primitive vector of Z[i] as an element of P"'^(C)). 

Suppose the circles have two distinct points in common. The new point 
can be viewed as a non-zero primitive vector x in the Z-span of u and iv. 
Then either x or ix is in the span of iv and u + v. We will show that neither 
is possible. 

If X is in the span of iv and u + v, then the two (non-lax) vertex lattices 
{u, iv) and (iv, u + v) have two independent vectors iv and x in common. 
Hence they have the sum in common, and the two circles contain at least 
three distinct points in common, and are therefore the same circle (up to 
orientation). Hence the two lax vertex lattices are equal up to orientation. 
But then v = niv for n G Z, which is not possible. 

So it must be that it is ix, and not x, that is in the span of iv and u + v. 
But then x is a common vector in (u, iv) and {—v,iu + iv). This implies the 
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existence of a Z-linear combination of u, v, iu, iv. By the Z[i]-independence 
of u and v, the coefficients of this combination must vanish; and this imphes 
a; = 0, which is contrary to our assumption. 

Therefore the two circles have exactly one point in common: they are 
tangent. □ 

Proposition 12.3. The circles of Si in bijection with the vertices of a 
chamber of the Apollonian city form a primitive strongly integral Descartes 
quadruple. Furthermore, every primitive Descartes quadruple appears in 
Si, up to a rigid motion. 

The proof relies on a lemma. 

Lemma 12.4. // A,B,C,D G Z have no common factor, A,C > 0, and 
AC = -\- , then there are u,v & Z[z] with gcd(u, f ) = 1 such that 

A = uu, B = , C = vv, D = — — — . (3) 

Conversely, if u,v G Z[i] are coprime, then A,B,C,D given by (3) are 
integers with no common factor. 

Proof. First, we show that there is a bijection 

F\Q{i)) = {[A,B, C, D]:AC = B^ + D^} C 

given by 



uv + uv _ uv — uv 



' 2 ' ' 2i 

In fact, a brief calculation verifies that the image of the map falls in X. An 
inverse map is given by 

[A, B, C, D] ^[A,B- Di] = [B + Di, C]. 

(At least one of these is always defined, and they are equal when both 
defined since AC = B'^ + D^.) 

Now assume that xu + yv = 1, for some x,y G Z[i]. Then 

xxuu + xyuv + xyuv + yyvv = 1, 

or equivalently, 

__ _ f uv + uv\ _ f uv — uv\ 

xxuu + [xy + xy) I ^ j ~'~ ^ ^"^^ " ^"^^ \ — 2^ — j ~'~ ^^^^ ~ 

Hence if u,v E Z[i] are coprime, then A,B,C,D have no common factor. 
Therefore if we start with A, B, C, D without a common factor, A,C > 0, 
AC = B"^ + D^, and we choose some M,f G Z[i] which are coprime such 
that [u, v] (-7- [A, B, C, D], then (3) must hold exactly. □ 
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Proof of Proposition 12. 3. Consider a chamber of the Apollonian city. Propo- 
sition 11.1 implies that the corresponding circles are strongly integral. From 
Proposition 12.2, the circles form a Descartes quadruple. Suppose 



^i=l^l, "^=(71, w = —u — V. 



where u and v form a basis for Z[i]^. 

Let ~ Im(C^^) give the curvature of the circle which is the image 

of M under , as in Proposition 11.1. Write ^ = a + bi 

and 6 = d — ci respectively. 

One checks that then the curvatures corresponding the standard ultra- 
basis vertices for u, v, w are given by 

bi = —{ad — be) = B{u, iv), (4) 

62 = + rf^ + {ad — be) = B{w, iu), (5) 

63 = (a + c)^ + {b + d)^ + {ad - be) = B{v, iw), (6) 

64 = + 6^ + {ad — be) = B{w — iv, iu — v). (7) 

It is apparent that any common factor of all the bi must divide the norms 
of a + 6i, c + di and (a + c) + (6 + d)i, which in turn implies that these 
three Gaussian integers have a common factor. Since m, iv form a basis for 
this cannot occur. Hence the bi are coprime, and the quadruple is 
primitive. 

Now consider curvatures 61, 62, ^a, &4 ^ ^ satisfying the Descartes relation 
(1), and having no common factor. Without loss of generality assume that 
fci + 64, &2 + ^4 > 0. Rearranging (1) gives 

^ (61 + 62 - &3 + bif + bl = {bi + bi){b2 + 64). 

Note that (1) also entails that &i + 62 — ^3 + &4 = (mod 2). By primitivity, 
the integers bi + 64, 62 + &4, &4 and (61 + &2 — &3 + &4)/2 have no common 
factor. Therefore we are in the case of Lemma 12.4, and there are some 
a,b,e,d E Z such that 

f h-\f &I+&2-&3 + &4 , , . .gX 

[a — bt)[—e — di) = h 04^, (8j 

N{a - bi) = 61 + 64, (9) 
N{-e-di) = b2 + b4, (10) 

such that a + bi and e + di are coprime. 

We may choose u to have second coordinate a + bi and v to have second 
coordinate d — ci, while the first coordinates may be chosen so that u and v 
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form a basis for This is possible because a + bi and d — ci are coprime. 

Let w = — M — V. Then, from (8)-(10), 

^ ^ = — (ac + bd), 64 = —(ad — be) 

61 + 64 = + 6^ 62 + &4 = + c?^- 

and so equations (4)-(7) hold, which demonstrates that we have found 
a chamber giving the necessary curvatures. Finally, any two Descartes 
quadruples with the same set of oriented curvatures must be related by a 
rigid motion. □ 

The following lemma is a variation on Lemma 12.4. 

Lemma 12.5. Suppose that A, B,C, D E Z[i] are such that AC = + 
and both A and C are divisible by 2. Then there are a,b,c,d G Z[i] such 
that 

A = 2cd, B = i{ad — be), C = 2ab, D = ad + be. 
Proof. We begin by showing there is a bijection 

P^(Q(z)) X P^(Q(i)) = {[A,B, C, D]:AC = B^ + D'^} c P^(Q(i)) 
given by 

[a, c] X [b, d] I— i- [2cd, i{ad — be)., 2ab, ad + be] 

A brief calculation verifies that the image of the map falls in X. An 
inverse map is given by 

[A, B, C, D] ^ [C, Bi-D]x [C, -Bi - D] = [-Bi -D,A]x [Bi -D,A]. 

(At least one of these is always defined, and they are equal when both 
defined since AC = B^ + D'^.) 

Now suppose one has A, B,C,D G Z[i] where A and C are divisible by 
2, and AC = B^ + D^. Let C = Bi ~ D and let C = -Bi - D. Since 
AC = CC, we can find Gaussian integers Ai,A2, Bi and B2 such that 

G = ACi, C' = A2C2, A = A^A2, C = CiC2. 

Then taking 

/^i 

(j (j 

a=7r=<-'2, c= — =Ai, b= — = Ci, d= — = A2. 

L'l Oi 02 ^2 

gives 

{2cd, i{ad - be), 2ab, ad + be) = {2A, 2B, 2C, 2D). 

Since 2 = (1 + i)(l — i) divides A and C, it divides cd and ab. We also 
have that ad and be are congruent modulo 2. Let ir = 1 + i. It must be, 
therefore, that vr divides a and c or else b and d. In the former case, we 
may replace a with a/n and c with c/vr; this now gives 

{2cd, i{ad - be), 2ab, ad + be) = ((1 - i)A, (1 - i)B, (1 - i)C, (1 - i)D). 
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The latter case is similar. Repeating the argument with vr' = 1 — i, we 
obtain the necessary a, b, c, d. □ 

Lemma 12.6. If Zi, Z2, z^, G are centre-curvatures of a primitive 
Descartes configuration in Si, then Zi = Z2 = z^ = Z4 (mod 2). 

Proof. By the last statement of Proposition 11.1, it suffices to show that 
whenever two distinct circles in Si are tangent, they are in bijection with 
elements of the same coset of PSL2(Z[i])/ PSL2(Z) in PGL2(Z[2])/ PSL2(Z). 
Without loss of generality, we can assume one of the circles is the real line. 
Then it suffices to show that any image of M under M e PGL2(Z[i]) which 
is tangent to M satisfies M G PSL2(Z[i]). Since the point of tangency 
corresponds to a primitive vector of Z[z]^ with coordinates in Z, we can 
also assume without loss of generality, by pre-composing with an element 
of PSL2(Z), that the point of tangency is fixed by M. Conjugating M by 
an element of PSL2(Z), we may assume the point of tangency is the origin. 
Then M has the form 

'a 

Since M has determinant a, it must be that a is a unit. By precomposing 

with an element of PSL2(Z), we can make /3 purely imaginary. Since 

is not on the real line, it must be that a is ±1, and we are done. □ 

A criterion for four circles to form a Descartes configuration has already 
been described in [J 2, Theorem 3.1] and [16, Theorem 2.2], and it subsumes 
the usual Descartes relation. We cite it here, but in another form (under a 
change of basis). Let the product of the centre and curvature of a circle be 
called the centre- curvature. 

Theorem 12.7 (Theorem 3.1 of [ ]). Let 
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Suppose we have four circles of oriented curvatures (61, 62, &3, ^4) (^nd centre- 
curvatures {zi, Z2, Z3, Z4) . Define the curvature-centre-conjugate matrix 



M 



Zl 
Za 



Z2 

zi 
zij 



If the circles form a Descartes configuration, then M has nonzero first col- 
umn and satisfies 

/O o\ 

M^S^RSM = 2 . (11) 
lo 2 0/ 
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Conversely, any solution M to (11) with nonzero first column, and conju- 
gate second and third columns is the curvature-centre-conjugate matrix of a 
unique ordered, oriented Descartes configuration. 

Corollary 12.8. If Zi, Z2, z^, G Z[i] are the centre- curvatures of four 
circles in a primitive Descartes configuration, then there are generically two 
possible choices for the collection of curvatures bi, 62, &3, &4- These solutions 
correspond to two Descartes quadruples of circles which are carried to one 
another by the map z 1— —2/z. 

Proof. By assumption, there is at least one primitive Descartes configura- 
tion having given centre-curvatures, hence at least one possible collection 
of bi. If the known, the bi must satisfy the three conditions corre- 

sponding to the entries in the first row of the matrix in (11). The matrix 

/I -1-1 -1\ 



S^RS = - 

4 



-1 1 -1 -1 
-1 -1 1 -1 

V-1 -1-1 1/ 



is full rank. The vectors {zi, Z2, z^, z^) and {z^,z^,z^,Z4) are independent 
since {Re{zi))f^^ and {lm{zi))f^i ^.re (otherwise the quadruple would con- 
sist of circles with collinear centres) . Therefore the two homogeneous linear 
conditions are independent. Hence there is a two-parameter family of solu- 
tions which must satisfy a homogeneous quadratic condition. The condition 
that the bi are integers with no common factor determines the two solutions 
uniquely. 

Furthermore, if one applies the map z i— )■ —2/z to a Descartes quadru- 
ple of circles, then one obtains another Descartes quadruple. By Propo- 
sition 11.1, the centre-curvatures do not change. If the two quadruples 
are distinct, the curvatures must change (since the centre-curvatures and 
curvatures together determine the quadruple). □ 

Proof of Theorem 12.1. Proposition 12.3 shows that all possible Descartes 
quadruples (and hence Apollonian circle packings) occur in Si, up to rigid 
motions. It also shows that to each chamber is associated a primitive 
Descartes quadruple. It remains to show the converse. 

Since by Proposition 11.1 all circles in Si are strongly integral, let 2:1, ^2, ^3, ^4 G 
represent the centre-curvatures of four circles of Si in Descartes con- 
figuration, and let 61, 62, ^s, &4 G Z represent the oriented curvatures, 
respectively. Let 

A I ^2 — ^3 + ^4 ^ , p, 

A = Zi+ Zi, B = , C = Z2 + Z/i, U = Zi. 

By Theorem 12.7, B^+D^ = AC and S is a Gaussian integer. Furthermore, 
the numbers R.e{zi) must all be congruent modulo 2, as must the numbers 
Im(zj), by Lemma 12.6. Hence A/2 and C/2 are Gaussian integers. 
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We are in the case of Lemma 12.5, and there are some a, /3, 7, 5 G Z[i] 
such that 

2a/3 = A, i{aS - 7/5) = B, 2-f5 = C, a5 + 7/5 = D. (12) 

Set 

^=(^)' "^=(3)' ^=-«-^- 

Let ^ ^ ~ ^^^^ " ^^^^ centre-curvature of the circle which 



V2 



is the image of M under 1 J ^ J , as in Proposition 11.1. One 

checks that 

Zi = 2af3 — a6 — 'yf3 = Z{v, iw), (13) 

Z2 = 275 — a5 — 7/3 = Z{w, iu), (14) 

= 275 + 2a(3 — 2*7/3 + 2ia6 — a6 — •-^(3 = Z{w — iv, iu — v), (15) 

Z4 = a6 + 'y(3 = Z{u,iv). (16) 
Also a consequence of Theorem 12.7 is that 

BB + DD- -AC - -AC = 2. 
2 2 

However, we have found expressions for A,B,C,D in terms of a,/3,7,5; 
plugging these in, we obtain 

BB + DD- ^AC - ^AC = 2NQ^iyQ{a6 - 7^), 

where iVQ(j)/Q is the Gaussian norm. This implies that a6 — 7/3 is a unit. 
Therefore, u and v form a Z[z]-basis of and the values zi, Z2, Z3, z^ are 
the values of Z at the corners of one chamber of the Apollonian city, i.e. 
we have found a Descartes quadruple with the required centre-curvatures. 

If we replace (a,/3,7,5) with (5, 7, /3,a) in (12), then we recover the 
values {C,B,A,D) instead of {A, B,C, D); this corresponds to the same 
set of centre-curvatures in a new order Z2, Zi, Z3, Z4^. Then we define instead 

f5\ A/3\ 

M= _ , f= ._ , W = —U — V. 

Then instead of (13)-(16) we have 

Z2 = 275 — q;5 — 7/3 = Z{v, iw), 
z\ = 2a(3 — a6 — 'yf3 = Z{w, iu), 

Z3 = 275 + 2af3 — 2*7/3 + 2ia6 — a6 — 'jf3 = Z{w — iv, iu — v), 
Z4^ = a6 + 'jf3 = Z{u, iv). 
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As before, u and v form a Z[2]-basis of Z[z]^. This represents a second cham- 
ber with the same centre-curvatures, but generically different curvatures. 
These two different chambers correspond to quadruples which are related 
by applying the map z —2/z on C. 

Corollary 12.8 implies that these two primitive Descartes quadruples cor- 
respond to the two possible quadruples with given centre-curvatures. There- 
fore, one of the two chambers we've found is exactly the quadruple we began 
with. 

An Apollonian circle packing is formed by completing all triples of circles, 
beginning with a specific originating quadruple; correspondingly, each of the 
completing circles is obtained as the vertex of an adjacent chamber. Hence 
all the circles in the packing correspond to vertices in the same palace as 
the originating chamber, and all the vertices in the palace give rise to circles 
in the Apollonian circle packing. □ 

13. Apollonian Circle Packings 
Consider the group 

r'= °VV ° 



ij \0 i 

This group acts by right multiplication on lax vertex lattice matrices. If 
we consider the palace containing the lax vertex lattice (1,0), (0, 1) corre- 
sponding to the real line, then T' takes a lax vertex lattice in this palace to 
all other lax vertex lattices in the same palace (by the proof of Proposition 
9.1). 

Hence the following is a corollary to Proposition 9.1 and Theorem 12.1. 
Corollary 13.1. The collection of images of the real line under the Mobius 
transformations ^^g^ ? ^ NT' for some N G PGL2(Z[z]) form a primitive 
integral Apollonian circle packing. 

Figure 7, which is obtained as all images of the real line under P', shows 
all the lax vertex lattices in the palace containing the lax vertex lattice 
(1,0), (0, 1). We obtain the horizontal Apollonian strip packing (it fills an 
infinite horizontal strip with this repeating pattern; the horizontal lines are 
circles of zero curvature). 

The curvatures of all the circles shown in this packing are even, which 
explains the convention of dilating^ these pictures by a factor of 2 in The- 
orem 8.1. Hence, we obtain the more usual strip packing (with an infinite 

string of unit circles) by taking all images of M under i ^ P'. 

^This is the lesser of several evils, each of which results in a proliferation of terrible 
twos. 




If we use instead a left coset A^F' of F', for some G PGL2(Z[z]), we 
obtain other Apollonian circle packings. Figure 8 is obtained upon using 

If we take together all the left cosets of F' for A^ G 

PSL2(Z[z]), drawing all the images of the real line, we obtain an Apollonian 
superpacking, as studied by Graham, Lagarias, Mallows, Wilks and Yan 
[13]. This picture, which those authors call the standard strongly integral 
super-packing, is shown in red in Figure 2. This repeating pattern contains 
one copy of every Apollonian circle packing in its unit square (except the 
strip packing). For details, see [13]. 

From the explicit bijection between lax circle lattices and certain equiv- 
alence classes of binary quadratic forms which we will see in Section 14, 
it will be apparent that a class number counting quadratic forms of given 
discriminant can be used to count the number of distinct occurrences of a 
given curvature in a region of the super-packing. This was already stud- 
ied in an earlier paper by Graham, Lagarias, Mallows, Wilks and Yan [11, 
Theorem 4.2]. 

The Apollonian superpacking has a dual, which consists of the images 
of the real line under all ^ NT' for N e (^^ PSL2(Z[z]), and 

is shown in blue in Figure 2. The superpacking together with its dual, 
consist of all images of the real line under i ^ PGL2(Z[i]). The 

resulting riot of circles can be considered a picture of the collection of full 
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-l+i 








Figure 8. The Apollonian (-1, 2, 2, 3) packing, obtained as 
all images of R under ( 1 ) ( j ° )r' . 

rank sublattices of (modulo lax equivalence), as described in the next 
section, or, simply, a picture of PGL2(Z[z])/PSL2(Z). 

14. The view from a vertex: Conway's sensual quadratic form 

The Apollonian city is not only formed in analogy to Conway's topo- 
graph, but it in fact includes infinitely many copies of the topograph. The 
topograph can be highlighted by choosing any one vertex of the Apollonian 
city and restricting attention only to those vertices at distance one from the 
chosen vertex. 

Theorem 14.1. Fix a vertex of the Apollonian city with vertex lattice L. 
Then the adjacent vertices may be put in bijection with the primitive ele- 
ments X E L in such a way that the following hold. 

(1) Two such vertices are connected by an edge of the Apollonian city if 
and only if their corresponding x,x' E L form a basis for L. 
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(2) Let H he the imaginary part of a Hermitian form. Then the vertices 
have H-values 

H{ix,x) ~ H{L). 

Furthermore, H{ix, x) is a binary quadratic form on the Z-module 
L of discriminant 

4 {H{u, vy - H{iu, u)H{iv, v)) , 

where u,iv is any ordered basis of L. 

Proof. Fix any one vertex of the Apollonian city, which we will call the 
mother vertex. Suppose it is labelled with the vertex lattice L generated 
by u and iv, and suppose that we have labelled this palace with one of the 
two possible global split labellings, as in the proof of Proposition 3.2. Now 
survey the motherland: those vertices at distance one from the mother, 
which we will call subjects, together with the edges which connect them. 

Fix a subject. Consider the split labelling on the edge joining the mother 
to the subject. The label closest to the mother is a vector x which is 
in the mother's vertex lattice. Thus we associate each subject with this 
corresponding x. 

Two vertices in the motherland are joined by an edge if and only if, 
together with the mother, they form a wall in the Apollonian city. This is 
the case if and only if together, they form a Z-basis for L (see the proof of 
Proposition 3.2). 

The function 

g : Z[z]^ — 7- Z, g{x) = H{ix,x) = h{x,x), 

satisfies the parallelogram law g{x + y) + g{x — y) = 2{g{x) + g{y)). Re- 
stricting g to L, we obtain a binary quadratic form. It takes values on the 
collection of edges connecting the mother with her subjects, so it can be 
viewed as a form on the motherland. The values it takes are exactly the 
if-labels of those edges in the Apollonian city. The discriminant of the form 
can be calculated as 

4 {H{u,vf - H{iu,u)H{iv,v)) . 

This association of a subject in the motherland to a primitive vector 
X E L implies that there exists some basis x, iy of L. Let z be such that 
X + y + z = 0. Then the subject has vertex lattice generated by z,ix. 
Property (H2) of H implies that 

H{z, ix) = g{x) — C 

where the constant C = H{x, iy) is the if-label of the mother. In other 
words, the if-labelling of the collection of vertices in the motherland is a 
translated quadratic form. 

Finally, note that the definition of the quadratic form is actually only 
dependent on the lax vertex lattice. □ 
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With reference to the language of the proof, suppose we illustrate just 
those edges of the Apollonian city which connect vertices in the motherland, 
and label the vertices by the vector of L associated to them. Then we obtain 
something very familiar (compare the first picture in Section 5): 

3n — iv 3u- 



2u — iv 



3u — 2iv 



2u — 3iv 




3u + 2iv 



2u + 3iv 



u — 2iv 



u + 2iv 



u — 3iv 



u + 3iv 



In other words, the motherland is a copy of Conway's topograph (in its 
'dual' triangular form). 



15. Circle lattices 

The curvature form 

5((«,/3),(7,5)) = Im(/35), 

of Proposition 11.1 depends only on the second coordinates of its arguments. 
Considering /3 = a + bi, 6 = c + di G Z[i] as vectors (a, 6), (c, d) e Z^, we can 
think of each vertex of the Apollonian city as being labelled by an oriented 
full-rank sublattice of generated by (a, 6) and (c, rf), instead of a vertex 
lattice generated by {a, /3) and (7, 6). This is a circle lattice, as in Definition 
4.3. The corresponding lax circle lattice is an equivalence class of two circle 
lattices: 

fa c\ f b d\ 
\b d) ' \-a -c) ' 

Proposition 15.1. Given a circle in Si, there is a unique lax circle lattice 
that may be associated to it. Furthermore, the circle has oriented curvature 
equal to the negative of the determinant of the circle lattice. 

A circle lattice determines the oriented curvatures of the Descartes quadru- 
ples and the Apollonian circle packing in Si containing any circle to which 
it is associated. 
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Proof. The first statement is a consequence of Theorem 8.1, since the circle 
lattice is determined by the vertex lattice. The circle lattices ^ and 



b d 

—a —c 



both have determinant 



ad — be = —B{{a, (3), (7, 6)), 



which is the negative of the curvature of the circle associated it. 

For the last statement, the bases of the circle lattice are the second co- 
ordinates of the bases u and iv of the associated vertex lattice. This, in 
turn, determines the curvatures of any Descartes quadruple associated to 
a chamber containing that vertex, by the proof of Proposition 12.3. In 
particular, if the second coordinates of u and v are a + bi and d — ci, respec- 
tively (corresponding to the form of the circle lattice above), then for the 
first standard ultrabasis chamber associated to superbasis u,v,w = —u — v, 
these are 

bi = —{ad — be) = B{u, iv), 

62 = + rf^ + {ad — be) = B{w, iu), 

63 = (a + c)^ + {b + d)^ + {ad - be) = B{v, iw), 

64 = ci^ + + {ad — be) = B{w — iv, iu — v). 

Since Descartes quadruples in Si are in bijection with chambers (Theorem 
12.1), this determines the curvatures of any Descartes quadruple containing 
a circle with that circle lattice. The curvatures of one Descartes quadruple 
determine the curvatures of the entire Apollonian circle packing of which it 
is a member. □ 

A given lax circle lattice may occur more than once in an Apollonian circle 
packing, but only if a given Descartes quadruple occurs repeatedly: the only 
example in which it may occur infinitely often is the strip packing, where, 
for example, each of the infinitely many circles of radius one necessarily has 
circle lattice Z^. All other Apollonian circle packings contain at most finitely 
many circles of curvature below a given bound (see [5, 11] for counting 
estimates) . 

The packing shown in Figure 8 has a symmetry, and therefore the same 
lax circle lattice often occurs twice or four times. In Figure 4, we show a 
circle lattice associated to each circle in that packing. 

It is nice to record the resulting lattice Descartes rule, which gives, for 
any three circle lattices which appear in an Apollonian circle packing, the 
circle lattices attached to the two circles completing the triple to a Descartes 
quadruple. 

Theorem 15.2. Consider three mutually tangent circles in an Apollonian 
circle packing. For some a, b, c, d, they have the following circle lattices 
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(represented by matrices whose column vectors form a basis): 

yb a J c —a — a J y—b + c a J 

The two completing circles have the following circle lattices: 

„ _ f —b — d c + a + d\ — a d — b + c \ 

^~\a + c d + b-cj' ^~\d-b -c + a + dj- 

Proof. By Theorem 12.1, the triple of circles is associated to a superbasis 
in the Apollonian city. Choose vectors u,v,w representing the superbasis 
such that u + V + w = 0, choose a global split labelling (as in the proof 
of Proposition 3.2), and hence a labelling by non-lax vertex lattices, where 
one vertex has vertex lattice generated by u, iv. For the first coordinates of 
u and iv, write a + bi and c + di respectively. The result follows. □ 

In practical application, to find the matrix forms for A, B, C above, so 
that one may apply this theorem, one may use the proofs of Proposition 
12.3 and Lemma 12.4. 

16. Quadratic forms in curvatures 

Besides a circle lattice, each circle also has an associated quadratic form, 
as described in Section 14. In this section we consider the quadratic form 
obtained in Theorem 14.1 by using the curvature form B as in the last 
section. 

Graham, Lagarias, Mallows, Wilks and Yan have already described a qua- 
dratic form associated to a root quadruple of an Apollonian circle packing 
in [11, Theorem 4.2]; their form lies in the equivalence class associated to 
the 'outer' circle of a packing in what follows. Bourgain and Fuchs have also 
described the same form in terms of a quadruple of curvatures [•>]. These 
are all essentially the same form. 

For a circle associated to a vertex lattice with basis u, iv having second 
coordinates /3 and 5 respectively, this quadratic form obtained from the 
curvature Hermitian form is 

g{Xu + Yiv) = iVQ(,)/Q(X/3 + Y6). 

Using the convention that P = a + bi and 6 = c + di, the discriminant of 
this form is related to the curvature of the circle: 

-AH{{a, P), (7, = -A{ad - bcf, 

A representative of the equivalence class is given by 

PPX^ + 2 Re{p5)XY + 

= {a^ + b'^)X^ + 2{ac + bd)XY + (c^ + d^)Y'^. 
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Figure 9. The Apollonian strip packing, with curvatures. 

The matrix form Q of this quadratic form is related to the matrix M of the 
circle lattice (whose column vectors are the first coordinates of u and iv). 
In particular, M^M = Q: 

fa b\ fa c\ _ f + 1)^ ac + bd\ 
\c dj \b dj ~ \ac + hd + d^ J ' 

Note that the resulting quadratic form is a primitive form: gcd(a^ + 6^, ac + 
hd, + d^) = 1 (otherwise a + hi and c + di are not coprime). 

The circles which are tangent to a fixed circle of lattice /3, 5 form the 
motherland, as in Section 14. They have curvatures 

A/'Q(i)/Q(X/3 + Y6) -ad + be. 

For example, consider the unit circle in the strip packing (Figure 7, drawn 

dilated by a factor of two). It has circle lattice ^ ^ . The circles tangent 

to it have values + — 1, for all primitive vectors {X,Y) G Z^. For 
example, one sees for all n G Z but — 1 only when n = 1 or n is the 
largest member of a pythagorean triple. These values are 

0, 1,4, 9, 12, 16, 24, 25, 28, 33, 36, 40, . . . 

See Figure 9. 

Proposition 16.1. If matrix 

represents a primitive positive definite quadratic form, whose discriminant 
is — 4n^ for some n G Z, then 

Q = M^M 
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for some integral matrix M of determinant ±n. Furthermore, there are 
exactly four solutions for M , of the form 

fa c\ f —a —c\ f h d\ f —h —d\ 
\h d)\-h -d)\-a -c) \a c)' 

Conversely, if 



is any integral lattice, then M'^M is the matrix of a primitive positive def- 
inite quadratic form with discriminant —4 (ad — bc)"^. 

Proof of Proposition 16.1. The existence of M amounts to solving the equa- 
tions 

A = a^ + h'^,B = ac + hd,C = c^ + D = ad - be, AC - B'^ = D^. 

Since the discriminant is negative, AC > 0. But since the form is positive 
definite, A,C > 0. From primitivity. A, B, C, D have no common factor. 
The result follows from Lemma 12.4. 

The converse is immediate. A form with negative discriminant but A,C > 
is positive definite. □ 

Hence a primitive quadratic form of the type of Proposition 16.1 arises 
from some circle lattice, and is part of a Descartes quadruple of some pack- 
ing. One can recover a Descartes quadruple simply be evaluating the qua- 
dratic form on a superbasis. 

Corollary 16.2. We have an explicitly computable bisection: 



The lattice Descartes rule of Theorem 15.2 gives rise to a corresponding 
quadratic form Descartes rule. 



Bestvina and Savin [2] study a complex formed of bases and superbases 
in Because a basis is contained in four superbases, each basis is 

associated to a square 2-cell, whose edges are the superbases. Square 2- 
cells are glued three along each edge, reflecting the fact that a superbasis 
contains three bases. The link of a vertex is the 1-skeleton of a cube, which 
describes which superbases share common bases. 

To compare Bestvina and Savin's complex to the Apollonian city studied 
here, consider the six faces of the cubical link. These correspond to vectors 





17. Appendix: Bestvina and Savin's topograph 
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u, V, w, u + iv, V + iw, w + iu, the edges of a tetrahedron in the Apollonian 
city. The eight vertices of the cube can be broken up into two collections 
of four, each of which forms the vertices of a tetrahedron: the collection of 
diagonals of each face which connect the four vertices gives the edges of the 
tetrahedron. 

In this way, each vertex of Bestvina and Savin's complex contains two 
dual tetrahedra. As mentioned briefly in Section 3, an ultrabasis tetrahe- 
dron has a natural dual, since the vectors labelling the edges adjacent to a 
vertex form a superbasis. This corresponds to the dual Descartes quadruple 
described by Coxeter ['■)]. 

There are other pieces of the vector/basis/superbasis/ultrabasis descrip- 
tion of P^(Q(2)) that remain to be studied; these complexes are only some of 
the possible complexes formed by the relationships between these objects. 

18. Final comments 

A future paper will study the relationships between the reduction of qua- 
dratic forms, continued fractions in the Gaussian integers, and reduction 
in Apollonian circle packings. Other natural questions include the study of 
other quadratic number fields (noteably the field with unit group of order 
six), as well as a wide range of questions concerning the arrangement of 
lattices into circle packings. 
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